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Turbulence measurements from H I absorption spectra 
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ABSTRACT 

We use the millennium Arecibo 21 cm absorption-line survey measurements to examine the 
issue of the non-thermal contribution to the observed Galactic H I line widths. If we assume 
a simple, constant pressure model for the H I in the Galaxy, we find that the non-thermal 
contribution to the line width, v n t scales as v^ t oc l a , for v n t larger than ~ 0.7 km s _1 . Here 
I is a derived length scale and a ~ 0.7 ± 0.1. This is consistent with what one would expect 
from a turbulent medium with a Kolmogorov scaling. Such a scaling is also predicted by 
theoretical models and numerical simulations of turbulence in a magnetized medium. For non- 
thermal line widths narrower than ~ 0.7 km s _1 , this scaling breaks down, and we find that 
the likely reason is ambiguities arising from Gaussian decomposition of intrinsically narrow, 
blended lines. We use the above estimate of the non-thermal contribution to the line width 
to determine corrected H I kinetic temperature. The new limits that we obtain imply that a 
significantly smaller (~ 40% as opposed to 60%) fraction of the atomic interstellar medium 
in our Galaxy is in the warm neutral medium phase. 
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1 INTRODUCTION 

The classical description of the Galactic atomic interstellar medium 
(ISM) is that it consists of the cold neutral medium (CNM) and 
the warm neutra l medium (WNM), in rough pressure balance wit h 
each other (e.g. lFiel j[l963 : iField et all 19691 ; IWolfire et alj|2003h . 
Detai led modeling of the energy balance in a multi-phase medium 
(e.g. IWolfireetal.ll 19951) shows that the pressure equilibrium can 
be maintained for H I in one of two stable ranges of kinetic temper- 
ature (T K ), viz. ~ 40K -200K for the CNM and ~ 5000K -8000 
K for the WNM. H I at intermediate temperatures is unstable and 
is expected to quickly migrate into one of the stable phases, unless 
energy is intermittently being injected into the medium. Recent ob- 
servations and simulations indicate that in our Galaxy a significant 
fraction of the atomic ISM is in the thermally unstable region (e.g. 
Vazquez-Semadeni et al]|200d : lHeilesll200ll ; iHennebelle & Audit! 
20071 - iHennebelle et al.ll2007l) . 

The classical method of determining the temperature of the 
atomic ISM is to compare the H 1 21 cm line in absorption towards a 
bright continuum source with the emission spectrum along a nearby 
line of sight. Assuming that the physical conditions are the same 
along both lines of sight, one can measure the sp in temperature (Ts) 
(or ex citation temperature) of the H I (see e.g. iKulkarni & H eiles 
1988, for details). While the H I spin temperature, strictly speak- 
ing, characterizes the population distribution between the two hy- 
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perfine levels of the hydrogen atom, it is often used as a proxy for 
the kinetic temperature of the gas. This is because, in high density 
regions, Ts is expected to be tightly coupled to the kinetic tempera- 
ture via collisions, while in low density regions, resonant scattering 
of Lyman-a photons is generally e xpected to c ouple the spin tem- 
perature to the kinetic temperature dFieldlll958l) . 

The 21 cm optical depth of the WNM is extremely low (typ- 
ically < 10~ 3 ) which makes it very difficult to measure the H I 
absorption from gas in the WNM phase. Consequently, emission- 
absorption studies usually provide only a lower limit to Ts. If the 
particle and Lyman-a number densities are low, Ts could in turn 
be significantly lower than Tk- On the other hand, one could use 
the 21 cm emission line width to determine an upper limit to the 
kinetic temperature. The line width is an upper limit to the temper- 
ature because in addition to thermal motions of the atoms, both bulk 
motion of the gas (e.g. differential rotation) as well as turbulence 
contribute to the observed line width. 

The presence of turbulence in the atomic ISM of our own 
Galaxy can be detected through, for e.g. the scale free nature of the 
power spectrum of the intensity fluctuat ions in H I 21 cm emission 
dCrovisier & Dicke yj [l983l ; lGreenlll993l) . In a turbulent medium, 
one would also expect the velocity dispersion to increase as a power 
of the length scale. Such a power law velocity width length scale 
scaling has been ob served in the ato mic ISM of the Large Magel- 
lanic Cloud (LMC) jKim et alj2007t) . To the best of our knowledge 
it has not been observed in the atomic ISM of our own Galaxy. 
In this paper we show that, assuming that the atomic ISM is in 
rough pressure equilibrium, the data from the millennium Arecibo 



2 N. Roy, L. Peedikakkandy and J. N. Chengalur 



- n. n ..^ n --.. n ..- 



log(N HI ) (1<r u crTO 

Figure 1. Nhi - Ts plot for Gaussian components from the millennium 
Arecibo 21 cm absorption-line survey data. The filled squares are compo- 
nents detected both in emission and absorption and the empty squares are 
components detected only in emission giving a lower limit of Tg . 



21 cm absorption-line survey jHeiles & Trolandl2003j|bT) is con- 
sistent with a velocity-length scale relation of the form a 2 oc l a . 
We also show that this scaling is, to zeroth order, consistent with 
that expected from turbulence in an medium with magnetic field of 
~ few /iGauss. 

Once one has an estimate of the turbulent velocity contribu- 
tion to the observed velocity width, one can correct for it, to derive 
a tighter limit to the kinetic temperature. We show that this correc- 
tion leads to a substantially smaller fraction of the gas being in the 
WNM phase than if one does not take turbulence into account. 



2 DATA, ANALYSIS AND RESULTS 

The data we use are taken from the millennium Arecibo 21 cm 
absorption-line survey and consist of emission and absorption spec- 
tra with a velocity resolution of ~ 0.4 km s -1 towards a total of 
79 background radio sources. The observational and analysis tech- 
niques are discussed in detail by iHeiles & Troland 1 2003al) and 
the astr ophysical implications are discussed in IHeiles & Troland 
< l2003bl) . A brief summary is that the absorption spectra were cor- 
rected for emission from gas in the telescope beam by interpolat- 
ing multiple off-source spectra, after which both the emission and 
absorption spectra were modeled as a collection of multiple Gaus- 
sian components. For each component, the spin temperature, upper 
limits on kinetic temperature, column densities and velocities were 
derived using these fits. There are several syst ematic uncertain- 
ties in s uch an analysis, discussed for example in Heil es & Troland 
(2003a), in particular those associated with estimating and subtract- 
ing the emission, and the assumption that each Gaussian compo- 
nent is a physically distinct entity. While a more robust measure- 
ment of the absor ption spectra can be done using interferometric 
observations (e.g. iKanekar e t al. 2003), here we work with the fit 
parameters provided as part of the survey. All the systematic uncer- 
tainties relevant to the Arecibo millennium absorption survey hence 
also apply to our results. 

The survey lists a total of 374 Gaussian components in the 
emission spectra towards the 79 continuum sources Out of these, 
205 components are also detected in H I absorption and have Ts 
measurements. For 21 of these components either the spin temper- 
ature had to be set to zero (by hand) in order to attain convergence 



of the fit jHeiles & Troland 2003a) or the upper limit of the kinetic 
temperature computed from the line width is less than the spin tem- 
perature (or the lower limit of the spin temperature). The derived 
parameters for these components are clearly unphysical, and we do 
not use them in our analysis. We are hence left with a total of 353 
Gaussian components consisting of 188 components detected both 
in emission and absorption and 165 components detected only in 
emission. In Figure (TJ is a scatter plot of the spin temperature Ts 
against the column density Nhi ; filled points are components that 
have been detected in both emission and absorption while empty 
points are components that have been detected only in emission. It 
is clear from the plot that the survey did not detect much gas with 
Ts < 10 K. 

For a homogeneous cloud at temperature T, the pressure 
P = nkT, with n = Nhi/£, where Nhi is the column density 
and I is the length of the cloud. Putting these together we have 
I =NmkT/P. For the CNM clouds detected in absorption, it is 
quite reasonable to assume that the kinetic temperature is the same 
as the spin temperature Ts . If we further assume that the pressure is 
roughly constant across clouds, then we have I ocNhi Ts. Though 
the density and temperature of neutral ISM vary over a few or- 
ders of magnitude, this assumption is justified because the pressure 
changes, in most of the cases, only by a factor of a few since the 
turbulence in the gas is at most transonic. Further, for these com- 
ponents, the non-thermal component of the line width is given by 
vf lt oc (TKmoi — Ts), where TKmoi is the measured line width 
of this component. In Figure ([2} we show a scatter plot of (TK mM 
— Ts) against NhiTs, as discussed above, to zeroth order this can 
be regarded as a plot of non thermal velocity against length scale. 
The solid line in the figure is a dual power law fit; at large length 
scales (log(NmTs) ^ 21.4±0.2) the power law index is 0.7±0.1, 
while at small length scales, the power law index is consistent with 
zero. The dotted line shows a fit which assumes that the measured 
TKmoi is larger than the true TKmoi by 60K; it provides a reason- 
able fit to the data over five orders of magnitude in Nhi Ts. The 
length scale corresponding to a pressure of 2000 cm~ 3 K, as well 
as the non thermal velocity width in km/s are also indicated in the 
figure. 

The correlation between cloud scale length and the non- 
thermal line width that we see at long length scales is consistent 
with a turbulence driven velocity line-width relation. The assump- 
tions we have made (viz. homogeneous clouds, constant pressure) 
are fairly naive, and are certainly not exactly valid in the current 
situation. The large scatter around the fit would be partly due to a 
break down of the assumptions (e.g. variations in the pressure) and 
partly due to measurement errors. If variations in pressure were 
the dominant contribution to the scatter, one would expect to see a 
systematic variation of the scatter plot for components in different 
Galactic latitude ranges. However, no such variation is seen. 

The ISM is known to have clumpy density and velocity 
structures an d is believed to be turbulen t at sc a les ranging from 
auto kpc toieter et ai] 1 19761 : iLarsorj 1 198 it iDeshpande et all 
2000). Incompressible hydrodynamic turb ulence leads to th e 
famous Kolmogorov scaling o% oc l 2 ^ 3 (Kolmo gorovl Il94l[) . 
similar to what we see at large length scales. However, the 
Galactic ISM cannot be modeled simply as an incompressible 
fluid. Recent theoretical studies and numerical simulations have 
investigated in details the turbulence of multi-phase medium 
jKovama & Inutsukdl2002l;lAudit & Hennebeiiell2005l;lGazol et all 
120051 ; Ivdzquez-Semadeni et al. 120061; iHe nnebelle et al.l 120070 . 
In some of these casesi Vdzquez-S emadeni e t alj 20061 ; 
IHennebelle etafl 120071) . synthetic H I spectra are computed 
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Figure 2. NjjiTg — TK ma i- Tg plot for Gaussian components from the 
millennium Arecibo 21 cm absorption-line survey data. Components de- 
tected both in emission and in absorption are only shown here. The solid 
line is a fit to the data using two power laws. The dotted line is the model 
with a 60 K overestimation of Txmax (see jj2]for details). 



to study the effect of turbulence. In general, these analytical 
and numerical works predict a Kolmogorov- like turbulence in 
two-phase neutral ISM. Henne belle et alj 1 120071) also report, based 
on simulation results, a power law scaling erj; oc I ' 8 consistent 
with our observation. 

Now, since fractional ionization couples the H I to the mag- 
netic field, the turbulence is expected to be magnetohydrody- 
nami c (MHD) in nature. Thou gh simple and ingenious models 
(e.g. iGoldreich & Sridharl 1 1995t) of incompressible MHD turbu- 
lence have been proposed, most of the insights into incompress- 
ible and compressible MHD turbulence again come from numer- 

i II -j ° 

ical simulations (Cho et al. 2002, and references therein). Models 

(like lGoldreich & Sridharll 19950 predict a Kolmogorov-like energy 
spectrum, E(k) oc k~ 5 ' 3 , for incompressible MHD turbulence 
and this is supported by bo th numerical simulations and observa- 
tions (see ICho etalJ J2002h for details). In case of compressible 
MHD turbul ence, Alfven modes are l east susceptible to damping 
mechanisms dMinter & Span gler 1997) and hence the energy trans- 
fer in Alfven waves is of major interest. Again, numerical sim- 
ulations show that the energy spectra of Alfven modes follow a 
Kolmogorov-like spectrum. 

In a situation where the bulk of the energy transfer is via 
Alfven waves, the non-thermal velocity dispersion Sv is related to 
the magnetic perturbation ampli tude SB and H I number density 
riH as Sv = SB/JAiYunnrnK dArons & Max|[l975l |Roshj||2007l) 
where fi = 1.4 is the effective mass of an H+He gas with cosmic 
abundance, ran is the mass of the hydrogen atom and it is usually 
assumed that SB ~ B. Using this relation, the magnetic field is 
found to be of the order of few fiG (column density weighted mean 
and median values are 1 1.7 and 10.2 /iG respectively) with no sig- 
nificant trend related to "cloud" size. We note that there are various 
uncertainties to the derived equipartition magnetic field. But our es- 
timate is broadly consistent with the observed magnetic field in the 
diffuse neutral ISM and matches, within a factor of 2, with the me- 
dian magnetic field estimated for a sub- sample of these compon ents 
using Zeeman splitting measurements dHeiles & Troland 2005). 

The break that is clearly seen in Figure l[2} requires some at- 
tention. This change in the power law index can not be explained 
just in terms of lower signal to noise on the physical quantities at 
low NhiTs end. However, as shown in the figure, the data are well 



fit by a model in which the line width is overestimated by about 
~ 60 K. There are three systematic effects that may contribute to 
the overestimation of the line width without much affecting Nhi 
and Ts: (i) the finite spectral resolution, (ii) blending of two or 
more narrow components and (iii) velocity (but not Ts) fluctua- 
tions in the gas within the Arecibo beam. The contribution from 
the first effect is quantified by estimating the width of a Gaus- 
sian signal after smoothing it to a spectral resolution of 0.4 km 
s" 1 and adding noise similar to that in the actual spectra. The ef- 
fect is found to be almost negligible because of the high spectral 
resolution. A similar numerical exercise with two Gaussian com- 
ponents was done to check the effect of blending of narrow com- 
ponents and ambiguities in Gaussian fitting. In this case, the effect 
is most significant when the blended lines are of comparable am- 
plitudes and have separations comparable to their widths. For ex- 
ample, blending of components with TKmai = 60 K (width of 
the Gaussian ~ 0.7 km s _1 ) with a separation of ~ 1.2 km s _1 
results in typically 20 - 30 K overestimation of TKinax. When the 
amplitudes of two Gaussian profiles are comparable, the line width 
is overestimated by upto ~ 60 K. The third possibility, that is, a 
fine scale structure in the v elocity (but not in the temperature) has 
been proposed earlier (e.g. lBrogan et alJl2005l ; lRov et al. ||2006) to 
expla i n the observed fine sc ale H I opacity fluctuations dDieter et al.l 
ll97dlCrovisier et al . 1985). Such velocity fluctuations within the 
Arecibo beam will also cause an overestimation of TKmoi. We, 
however, note that the scale length (inferred from Nhi and Ts) of 
the components below the break is very small. Although the ex- 
istence of tiny "c louds" is supported by observations and numer- 
ical simulations teraun & Kanekarl 120051; Stanimirovic & Heiles 
120051: 



iNagashima et al 



| 2006l ; 



Vdzquez-Semadeni et al. 



l200d 

iHennebelle & Auditl 120071 . e.g.), their origin and physical proper- 
ties are still unknown. The evaporation timescale for these clouds 
are ~ 1 Myr. These structures can survive if either the ambi- 
ent pressure around the clouds is much higher than the standard 
ISM pressure or they are formed continuously with a compara- 
ble timescale. While we have presented plausible arguments for 
the break that we see not corresponding to a physical phenom- 
ena, the lack of detailed understanding of these tiny H I structures 
means that we can not rule out the possibility of some physical 
phenomenon being responsible for the break. 

2.1 A new indicator of the temperature 

For a multi-phase medium if the turbulent velocity dispersion scal- 
ing is similar for coexisting phases, then this scaling relation can 
be exploited to get a handle on the physical temperature of the gas 
that is detected only in H I emission but not in absorption. Since 
one only has a lower limit on Ts for these components, they lie, 
as expected, systematically on the top left side of the fit to the 
components detected in both emission and absorption (Figure $3$). 
For these components we define a proxy temperature Tl that will 
restore the component back to this power law correlation. Given 
the measured Nhi and TKmoi from the emission spectra, one can 
uniquely compute this proxy temperature. Since Tl corresponds to 
the velocity width after correction for the turbulent velocity, it is 
a better estimate of the actual physical temperature of the cloud 
than that of Txmoi • Note that since most of the components in Fig- 
ure ([3} line beyond the break in the fitted function the derived Tl is 
independent of the whether the break arises due to some underlying 
physical reason. 

For all except 2 of the 165 components detected only in emis- 
sion, Tl was calculated as described above. For two components, 
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Figure 3. NjjiTg — TK ma i- Tg plot for Gaussian components from the 
millennium Arecibo 21 cm absorption-line survey data. Components de- 
tected only in emission and not in absorption are shown here and Tg is the 
lower limit for these components. The solid line is derived from the com- 
ponents detected both in emission and in absorption. 
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Figure 5. Njji fraction, for components detected both in emission and in 
absorption (hatched histogram) and components detected only in emission 
and not in absorption (empty histogram), in different temperature range us- 
ing Tl and Tn max as the proxy for the actual physical temperature. 
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Figure 4. Histogram of derived temperatures Tl and Tk ■ The solid lines 
only include components detected both in emission and in absorption and 
the dotted lines include all the components. The top two panels give the 
number of Gaussian components and the bottom panels give the H I column 
density. The left two panels are histogram of Tl and the right two are of 
T^max. 



no meaningful solution for Tl could be found, and they are hence 
not included in the further analysis. For the components detected 
both in emission and in absorption, Ts is taken to be same as Tl. 
With this, we have Tl and Tumax for a total of 351 Gaussian 
components. Figure 10 shows the histogram of Tl and TKmoi in 
terms of both number of "clouds" and H I column density. The 
top two panels give the number of Gaussian components and the 
bottom panels give the H I column density in different tempera- 
ture bin. From the histograms, it is evident that our results qualita- 
tively confirm the earlier detection of a significant fraction of gas in 
the thermally unstable region. Quantitatively, however, a significant 
fraction of the gas with high Tumax after correction for turbulent 
broadening corresponds to gas in the stable phase. This quantitative 
difference is illustrated in Figure ((5} which shows the Nhi fraction 
for both the population (components detected both in emission and 
in absorption and components detected only in emission) in dif- 
ferent temperature range using Tl instead of Tktmi as the proxy 
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Figure 6. Histogram of derived temperatures for Tl < 500 K. The solid 
lines only include components detected both in emission and in absorption 
and the dotted lines include all the components. 



for actu al physical temperature. Further, while iHeiles & Trolandl 
d2003bl) find that 60% of all H I is in the WNM phase, we find 
that only 40% of the gas has temperature higher than 500 K. 

A closer examination of the components with Tl < 500 
shows a clear peak near Tl ~ 50 K in the number distribution and 
that the major fraction of the gas is below Tl ~ 100 K as shown 
in Figure l[6j, Figure Q shows the histogram of scale length L ~ 
Nhi Ts. This clearly shows a bi-modal statistical distribution of the 
"cloud" size for the neutral ISM. As expected, the dominant contri- 
bution to the the peak at lower L is from the cold components and 
to the peak at higher L is mostly from the warm components. If N 
is number of clouds along the lines of sight, R is typical size of the 
clouds and n is H I number density, then N(H l)x oc NxRxnx 
where X stands for CNM or WNM. From the observed N(H I) and 
N for all the components used for this analysis we find, using this 
relation, -Rwnm/-Rcnm ~ 110 for typical ticnm/^wnm ~ 100. 
This is consistent with the ratio of the length scale corresponding 
to two peaks in the bi-modal distribution. 
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Figure 7. Histogram of scale length L ~ Nhi Ts in terms of number of 
components (solid line) and total H I column density (dotted line). Two 
distinct peaks are evident in both cases. 



3 CONCLUSIONS 

In this work we present a new phenomenology based technique to 
address the issue of non-thermal line width and the temperature of 
the diffuse neutral hydrogen of our Galaxy assuming a rough pres- 
sure equilibrium between different phases of the ISM. A possible 
connection between the observed Kolmogorov-like scaling of the 
non-thermal velocity dispersion in the Galactic H I and the turbu- 
lence of the interstellar medium is discussed. This scaling relation 
is used to re-examine the issue of the temperature of the Galactic 
ISM with the help of the millennium Arecibo 21 cm absorption- 
line survey measurements. The distribution of the derived temper- 
ature is found to be significantly different from the distribution of 
the upper limits of the kinetic temperature. A considerable fraction 
(~ 29%) of the gas is found to be in the thermally unstable phase, 
qualitatively confirming earlier results. However, about 60% of all 
the neutral diffuse gas, a much higher fraction than that of reported 
earlier, has temperature below 500 K. The CNM temperature dis- 
tribution shows a clear peak near T ~ 50 K and the cloud size for 
the neutral ISM shows a bi-modal statistical distribution. Derived 
magnetic field from the non-thermal velocity dispersion matches, 
within a factor of 2, with the magnetic field value estimated from 
the Zeeman splitting measurements. The Kolmogorov-like scaling 
is consistent with the existing theoretical prediction, numerical sim- 
ulations and earlier observational results. 
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